Diurnal change of retinal vessel density related to hemodynamic variation in treatment-naïve low-teens normal-tension glaucoma

This study undertook to investigate the diurnal variation of optical coherence tomography angiography (OCTA) -derived retinal vessel density (RVD) in glaucoma patients with low baseline intraocular pressure (IOP). A prospective evaluation was performed on low-teens normal-tension glaucoma (low-teens NTG) patients with pre-treatment IOP < 15 mmHg and 32 healthy subjects. Superficial peripapillary and macular RVD by OCTA, IOP, and systemic blood pressure (BP) were all measured four times per day (from 9:00 a.m. to 6:00 p.m.). In the low-teens NTG group, the magnitude of diurnal changes in peripapillary RVD and macular RVD were greater than those in the healthy group. Diurnal variations of diastolic BP (DBP) and mean ocular perfusion pressure (MOPP) also were greater in the low-teens NTG group. As for the patterns of diurnal RVD change, the inferior and temporal sections of macular RVD showed significant differences between the two groups. Diurnal changes of RVD and MOPP and were greater than those in healthy eyes. The macular RVD and MOPP showed different diurnal patterns between the two groups. From these findings, OCTA-derived RVD variation could be related to hemodynamic variability in low-teens NTG.

Intraocular pressure (IOP) being, at present, the most significant risk factor and sole reliable therapeutic target for the currently available treatments all aim to reduce IOP by pharmacologic, laser, or surgical intervention 1 . However, certain patients suffer from glaucoma even though they have relatively low IOP. In such cases, involvement of risk factors other than elevated IOP might be central 2 . The vascular theory, meanwhile, ascribes the pathogenesis of glaucomatous optic neuropathy (GON) to reduced perfusion pressure and faulty vascular autoregulation 2 . Compromised optic nerve perfusion likely plays a role in GON, and it may have a more central role in the glaucoma setting despite low or normal IOP 3 .
A recent study conducted by the present authors showed that fluctuation of diastolic blood pressure (DBP) and diurnal IOP variations were significantly associated with greater probability of disease progression in lowteens NTG eyes 4 . Based on these results, we could speculate that hemodynamic variability in low-teens NTG may play an important role in disease course. Meanwhile, advances in imaging devices such as optical coherence tomography (OCT) and OCT angiography (OCTA) have allowed for detailed visualization of superficial and deep vasculature as well as optic nerve head (ONH) structural characteristics 5,6 . In an earlier study, the present authors, having analyzed diurnal variation in choroidal thickness using OCT, reported that diurnal variation of choroidal thickness representing the deep vasculature was significant in patients with primary open-angle glaucoma (POAG) 7 . Subsequently, others noted that in open-angle glaucoma (OAG) eyes, diurnal change of superficial retinal vessel density (RVD) was significantly greater when compared with healthy eyes and that such change was correlated with mean ocular perfusion pressure (MOPP) change 8 . All of these findings suggest that diurnal RVD changes might reflect clinically important hemodynamic variation thought to be peasant in OAG 9 .

Scientific Reports
| (2023) 13:10613 | https://doi.org/10.1038/s41598-023-37214-w www.nature.com/scientificreports/ Those studies, however, had enrolled patients who were taking ocular hypotensive medication. Also, those participants usually were POAG patients with higher baseline IOP. In other words, our previous studies may not have fully represented actual hemodynamic variability in normal-tension glaucoma (NTG) patients, especially those not on ocular hypotension medication.
The present study therefore undertook to investigate diurnal variation of RVD and the related ocular and hemodynamic factors including IOP, blood pressure (BP), and MOPP in patients diagnosed with low-teens NTG compared with healthy subjects. The aim was to elucidate the pattern and extent of hemodynamic instability observable among patients with low-teens NTG.

Results
The baseline characteristics of the study population are presented in Table 1. There were no significant intergroup differences in age, sex, body-mass index (BMI), diabetes mellitus, hypertension, spherical equivalent, axial length, central corneal thickness, or baseline IOP (all P > 0.05). There also were no differences in number or type of antihypertensive treatment. As expected, low-teens NTG eyes had a higher vertical cup-to-disc ratio, a thinner average retinal nerve fiber layer (RNFL) thickness and worse visual field results than did the healthy eyes (all P < 0.001).
Assessment of repeatability of RVD. The repeatability of RVD obtained by SS-OCTA was excellent for the low-teens NTG patients and healthy subjects ( Table 2). The intraclass correlation coefficients (ICCs) for intra-visit repeatability of peripapillary RVD showed excellent reliability (0.835-0.948), as did those for macular RVD (0.863-0.948). The ICCs for inter-visit repeatability of peripapillary RVD ranged from 0.755 to 0.911, and for macular RVD, from 0.806 to 0.898.
Intra-group assessment of diurnal variations. The RVDs of the peripapillary (ONH, superior, inferior, temporal, nasal) and macular (fovea, superior, inferior, temporal, nasal) areas were assessed in detail ( Table 3). The data also was schematized as a scatter plot ( Supplementary Fig. 1). In the healthy group, only RVDs of the inferior peripapillary region showed significant diurnal variation (P = 0.008), whereas the low-teens NTG group showed sizeable variation in the superior peripapillary, inferior macular, and temporal macular areas (P = 0.018, Table 1. Demographics and clinical characteristics of healthy subjects and low-teens NTG. Values are presented as the mean ± standard deviation or number (%). *Student t-test, † Chi-square test. ACEI = angiotensin-converting enzyme inhibitor; BRB = beta receptor blocker; CCB = calcium channel blocker; cpRNFL = circumpapillary retinal nerve fiber layer.

Magnitudes of diurnal variations.
To determine the degree of diurnal variation, the maximal change (maximum-minimum values among 4 measurements) of RVD and hemodynamic parameters was analyzed ( Table 4). The magnitudes of diurnal variation of average peripapillary and macular RVD variations in low-teens NTG (3.85 ± 2.83 and 4.92 ± 3.83%) were more variable than in the healthy group (2.64 ± 1.55 and Table 2. Intra-and inter-visit repeatability of RVD in healthy subjects and low teens-NTG. *As analyzed by repeated-measures ANOVA between Healthy and POAG groups, respectively. ICC = intraclass correlation coefficient; CI = coefficient interval; RVD = retinal vessel density; ONH = optic nerve head.  .68 ± 9.67, P = 0.008), and nasal macular (5.09 ± 3.93 vs 8.38 ± 6.49, P = 0.007) areas in the low-teens NTG eyes were larger than in the healthy eyes. Similarly, among the hemodynamic variables, the maximal changes of DBP (12.65 ± 7.17vs 16.90 ± 7.35, P = 0.013) and MOPP (7.14 ± 4.43 vs 9.57 ± 5.27, P = 0.025) in low-teens NTG were greater than those in the healthy group.
Patterns of diurnal variations. The patterns of diurnal changes of peripapillary RVD ( Fig. 1), macular RVD ( Fig. 2), and hemodynamic parameters ( Fig. 3) between the two groups were examined, and the intra-and inter-group variabilities were assessed. The diurnal RVD of the peripapillary areas did not show any significant pattern difference between the two groups (all P > 0.05, Fig. 1). However, a comparative analysis of the diurnal macular RVD changes between the two groups found differences in some sectors. Specifically, the inferior and temporal macular RVD showed different patterns in the inter-group analysis (P = 0.029 and 0.036, Fig. 2). As for the hemodynamic parameters, the low-teens NTG patients and healthy subjects showed different patterns of diurnal variation in DBP and MOPP (P = 0.001 and 0.010, Fig. 3).
Post hoc analysis of pattern changes. For the variables showing inter-group differences, a post hoc analysis was performed to confirm which data-acquisition session(s) were most contributory ( Table 5). The inferior macular RVD of the healthy group was considerably higher than that of the low-teens NTG at 9 a.m.(P = 0.005), and temporal macular RVD was higher at 6 p.m.(P = 0.023). As for the hemodynamic parameters, the DBP and MOPP of the low-teens NTG group had trough values at 6 p.m., whereas the healthy group showed its peak DBP and MOPP values in that session (P = < 0.001 and 0.012).

Discussion
This study was a comparative analysis of the magnitude and pattern of diurnal RVD variation between low-teens NTG patients and healthy eyes. We found that the magnitude of diurnal RVD variation was more significant in low-teens NTG. As for the pattern of diurnal variation, macular RVD showed different patterns between the two groups. Characteristically, both DBP and MOPP also showed higher diurnal variation in the low-teens NTG group, with a descending pattern in the evening.
The present study has some strengths in its design. First, we studied exclusively low-teens NTG patients, who would be expected to be more vulnerable to hemodynamic variability compared with patients experiencing glaucoma with higher IOPs. Second, only treatment-naïve patients were studied; thus, we could exclude the effect of IOP-lowering medication on diurnal fluctuation. Third, IOP, RVD and systemic BP were simultaneously measured through a prospective analysis.
As mentioned in the introduction, the authors have already reported on the diurnal variation of RVD, IOP, and BP in POAG patients 8 . An earlier study of ours demonstrated that diurnal IOP variation in low-teens NTG is an independent risk factor for glaucoma progression 4 . According to multivariate regression analyses, the Table 4. Magnitude of diurnal variation of RVD and hemodynamic parameters between healthy subjects and low-teens NTG. *Independent t-test, bolded values represent significance, P < 0.05. RVD = retinal vessel density; ONH = optic nerve head; SBP = systolic blood pressure; DBP = diastolic blood pressure; IOP = intraocular pressure; MOPP = mean ocular perfusion pressure. www.nature.com/scientificreports/ significant risk factors for progression in low-teens NTG were fluctuation of DBP (HR = 1.058), diurnal IOP fluctuation (HR = 1.609), and occurrence of disc hemorrhage (HR = 3.661). According to these findings, diurnal variation of BP, in addition to IOP, is an important parameter to evaluate in cases of low-teens NTG. However, these patients with POAG mostly had high baseline IOP, and thus, mechanical damage from high IOP rather than hemodynamic factors might have been the dominant mechanism of GON. Also, because the participants had been using glaucoma medication, our previous studies may not have fully represented the actual hemodynamic variability. In this respect, this study was conducted to supplement the previous research, and the diurnal variation of RVD and MOPP before treatment were analyzed in low-teens NTG.
Notably, low ocular perfusion pressure (OPP) can occur as a result of low BP, even with normal IOP (< 21 mmHg). Also, in terms of OPP fluctuation, OPP sometimes is more sensitive to BP changes than to IOP changes, because BP varies in the range of 40 to 60 mmHg while IOP varies in the range of 5 to 8 mmHg 10 . In this regard, several epidemiological studies have identified low systolic BP (SBP) and DBP as an important risk factor for development and progression of glaucoma [11][12][13][14] . Some studies have also identified nocturnal BP dips as a risk factor for glaucoma progression 15,16 .
In particular, the DBP behavior of low-teens NTG in this study was distinct: it appeared at the same level as the healthy controls in the morning and afternoon (from 9 a.m. to 3 p.m.), but became significantly lower during the evening (6 p.m.), which coincided with the time during which the MOPP was lower. This would explain why MOPP was consistently lower in low-teens NTG in the evening. The association between DBP reduction and lower MOPP in low-teens NTG may be attributable to dysregulation of the autonomic nervous system 17 . A longitudinal study of chronic hypotension found that low BP cases were associated with low BMI, heart rate and body temperature, reflecting the lowered sympathetic nervous system in hypotension 18 .
In this study, we focused on not only MOPP level but also MOPP variability. The degrees of diurnal change of macular RVD accompanied by DBP and MOPP change were greater in the low-teens NTG group than in the healthy subjects. Previous studies have identified larger variability of MOPP as a risk factor for glaucoma progression 19,20 . Similarly, the present study observed significant variations of systemic hemodynamics and ocular hemodynamic changes simultaneously at the level of the retinal blood vessels in the low-teens NTG setting.
Under normal physiological conditions, low perfusion pressure is compensated by low resistivity to flow 21 . One possible explanation for this phenomenon is the effect of compensation mechanisms, which is to say, autoregulation. In healthy individuals, significant fluctuations in retinal circulation did not occur, as autoregulation to compensate for variation of MOPP works properly. Another OCTA study of ours on the systemic BP and IOP of 132 healthy participants (264 eyes) under 45 years of age demonstrated that peripapillary RVD was not correlated with MOPP 22 . Considering the young age and disease-free state of our study population, www.nature.com/scientificreports/ autoregulation most likely was at work. On the other hand, it could also be deduced that this compensation mechanism does not work properly in low-teens NTG groups 19,20 . If OPP is not constantly maintained and fluctuates markedly, oxidative stress on the ONH can result. OPP fluctuation can cause ischemia-reperfusion damage to retinal ganglion cells (RGCs) by oxygen-free radicals and nitric oxide species 23 . Additionally, OPP that fluctuates below the lower autoregulatory limit is more susceptible to RGC damage and ONH ischemia 24 . Verticchio Vercellin AC et al. 25 evaluated diurnal variations in peripapillary RVD as assessed by OCTA in healthy subjects as well as ocular hypertension and OAG patients. The results demonstrated that the variations in the OCTA-derived parameters were relatively small among the three groups. This discrepancy from our results appears to have been due to a difference of study design. Our study was conducted on low-teens NTG in the treatment-naïve state, whereas Verticchio Vercellin AC's study had evaluated diurnal variation in glaucoma patients being administered ocular hypotensive medication. According to previous reports, ocular hypotensive medication reduces diurnal fluctuation 26,27 . Therefore, the actual severity of diurnal RVD variation might have been underestimated in that earlier study.
In this study, significant diurnal variation of RVD was observed in the superior peripapillary sector, inferior macular sector, and temporal macular sector of low-teens NTG patients. There are some possible explanations for the sectoral differences in RVD variability seen in NTG eyes. The RNFL thickness in the superior and inferior peripapillary areas was relatively thicker than in the nasal and temporal areas, and the RVD was higher 28 . In addition, circumpapillary retinal vessel RVD was correlated with RNFLthickness 29 . That is, there is a possibility that the diurnal variation was measured larger in that area because the initial value was thicker. The next hypothesis that can be considered is that initial RNFL damage occurs in the superotemporal or inferotermpoal area, and especially the interotemporal area connecting to the macula is known as a macular-vulnerable zone. This zone is part of the high-density axon region known to be vulnerable to glaucomatous damage 30 . A follow-up study that topographically analyzes whether these glaucoma-damage-prone areas are correlated with higher variation of RVD seems necessary.
Interestingly, the RVD in the inferior region of the macula was lowest in the morning (9:00AM), not in the evening (6:00PM). Although a clear mechanism could not be delineated, the authors considered two possible answers. The standard deviation of the inferior macular RVD value at 9:00 am was relatively larger (± 7.56). And since this study had a relatively small sample size, we assumed that some outliers had made for low RVD at these points. Next, the 6 PM test time specified in this study may not be sufficient to reflect the actual rhythm of RVD change. The RVD value might decrease further after 6 PM, but in this study, the OCTA scans were conducted in the early evening (6 PM). www.nature.com/scientificreports/ Several limitations need to be considered when interpreting the results of this study. First, RVD and hemodynamics were not assessed at midnight. Uniformly across the entire diurnal cycle, notably, it was not practical to measure patients at the midnight time point. Thus, the significant pattern of RVD variation that we observed might not reflect the whole diurnal variation. Also, changes in IOP, BP, and RVD according to postural change were not considered in this study. In the near future, additional assessments are needed to analyze RVD based on 24-h monitoring of BP and IOP according to postural change. Second, another possible limitation that may have influenced diagnostic validity is 'white coat syndrome,' 31 which refers to stress on study participants in cases where their BP measurements are obtained by clinicians, which readings might be rendered inaccurate as a result. Third, MOPP, as an estimate of OPP, might differ from actual OPP at the ONH, for which reason brachialarterial pressure has been used as a surrogate, albeit not an ideal one, for ophthalmic arterial pressure 22 . Fourth, and finally, diurnal variations might be more or less pronounced, or have a different rhythm, in NTG. It might simply be that the rhythm could not be aligned with the specific times tested between the groups. Additional measurements for all subjects and on consecutive days would be helpful in clarifying this point.
In conclusion, in low-teens NTG eyes, diurnal changes of RVD, DBP and MOPP were greater than in the healthy eyes, and the patterns of diurnal variation significantly differed. Our findings suggest that diurnal RVD variation might reflect the hemodynamic variation characteristic of low-teens NTG. It is assumed that the diurnal variation of systemic BP and MOPP rather than of IOP will act as an important factor in diurnal RVD variation. Therefore, in addition to assessment of IOP fluctuation in low-teens NTG patients, diurnal variation of BP and MOPP should be considered when evaluating progressive low-teens NTG patients. In the future, complementary studies should be conducted to determine whether diurnal variation of RVD shows a significant association with functional and structural data and whether it can be a risk factor for glaucoma progression in NTG patients.

Methods
Ethics. This study was approved by the Institutional Review Board of Hallym University Sacred Heart Hospital, and all of the protocols were in accordance with the tenets of the Declaration of Helsinki. Each subject was required to sign an informed consent statement before being enrolled in the study and prior to any study measurements.

Study subjects.
This was a prospective cohort study of low-teens NTG and age-matched healthy subjects.
All of the subjects were prospectively enrolled at the Glaucoma Clinic of Hallym Sacred Heart Hospital between November 2019 and July 2021. Low-teens IOP was defined as mean untreated office IOP < 15 mmHg in diurnal measurements (at 3-h intervals from 9:00 a.m. to 6:00 p.m., for a total of 4 measurements). Low-teens NTG patients who met the following inclusion criteria were initially enrolled: (1) treatment-naïve status (diurnal variation analysis prior to antiglaucoma medication administration), and (2) normal open-angle on gonioscopy. All of the participants had to meet the following further criteria for inclusion: age > 18 years, best-corrected visual acuity (BCVA) ≥ 20/40, and spherical equivalent < -6.0 diopters. Patients with a history of ocular surgery except uncomplicated cataract surgery, or any history of retinal pathology or non-GON were excluded. Participants also were excluded in cases of diagnosis of neurodegenerative disease or cerebrovascular accident history. After evaluating the enrollment of low-teens NTG group, we then consecutively recruited healthy volunteers from a comprehensive ophthalmology clinic in order to match the age distribution of the glaucoma patients.
Participants with diabetes mellitus and systemic hypertension were included unless they had been diagnosed with diabetic or hypertensive retinopathy. Antihypertensive medications were recorded on the basis of the number of medications taken and the type of medication (angiotensin-converting enzyme inhibitors[ACEIs], angiotensin II receptor blockers, calcium channel blockers [CCBs],diuretics or beta receptor blockers [BRBs]) 32 .
All of the participants underwent complete ophthalmologic examinations, including baseline IOP measurement by Goldman applanation tonometry, BCVA, refractive value with an autorefractor (HRK-8000A, Huvitz, Gunpo, Korea), slit-lamp examination, gonioscopy, and dilated fundus examination. After pupil dilation, all of the participants were imaged by stereo optic disc photography and red-free (RNFL photography (Vx-10; Kowa Optimed Inc., Tokyo, Japan). They also undertook standard automated perimetry using the Swedish interactive threshold algorithm with the 24-2 standard program (Humphrey Field Analyzer II; Carl Zeiss Meditec Inc., Dublin, CA, USA).
Glaucomatous eyes were defined as those showing glaucomatous optic disc neuropathy (e.g., notching, neuroretinal rim thinning, and/or RNFL defects) and corresponding glaucomatous visual field defects, as confirmed by at least two consecutive visual field examinations. The glaucomatous visual field defects were defined as a cluster of ≥ 3 points with P < 0.05 on the pattern deviation map in at least one hemifield, including ≥ 1 point with P < 0.01; a pattern standard deviation (PSD) of P < 0.05; or glaucoma hemifield test result outside the normal limits. Glaucoma clinical severity was assessed for each eye based on mean deviation (MD) thresholds from the Table 5. Post hoc analysis of macular RVD, DBP, and MOPP in healthy group and low teens-NTG. *Independent t-test; bolded values represent significance, P < 0.05. RVD = retinal vessel density; DBP = diastolic blood pressure; MOPP = mean ocular perfusion pressure. www.nature.com/scientificreports/ modified Hodapp-Anderson-Parish classification scheme, as follows: early (MD ≥ − 6.00 dB), moderate (− 6.01 to − 12.00 dB), advanced (− 12.01 to − 20.00 dB), and severe (− 20.00 dB) glaucoma stages, respectively 33 .

Assessment of diurnal variation of RVD by swept-source optical coherence tomography angiography (SS-OCTA).
The OCTA scans were performed at four different time points over a single day or on two consecutive days, at 3-h intervals: 9:00 a.m., 12:00, 3:00, and 6:00 p.m. All of the scans were performed by a single operator (Supplementary Fig. 2).
In this study, we applied a mode of swept-source OCTA (SS-OCTA; DRI OCT Triton, Topcon, Tokyo, Japan) to visualize the retinal vasculature of the peripapillary area and macula. The scans were taken from a 4.5 × 4.5 mm cube, each consisting of 320 clusters of four repeated B-scans centered on the ONH and fovea. Slabs of retinal capillary plexus were segmented by the built-in software (IMAGEnet6, v1.23.15008, Basic License 10). Only the superficial capillary plexus (SCP) was adopted for analysis. The SCP was delineated from 2.6 μm below internal limiting membrane to 15.6 μm below the junction between inner plexiform and inner nuclear layers. The DRI OCT Triton characterizes vascular information quantitatively by vessel density (VD; %) as well as qualitatively at various user-defined retinal layers on a vessel density (VD) map. RVD was automatically calculated as the proportion of the area occupied by flowing blood vessels defined as pixels having above-threshold-level decorrelation values acquired by the OCTA RAT M (OCTA Ratio Analysis) algorithm 34 . Peripapillary and macular RVD were measured at the optic disc and fovea and at each of the adjacent four sites employing a modified version of the Early Treatment Diabetic Retinopathy Study (ETDRS) grid ( Supplementary Fig. 3).
Assessment of oculo-and hemodynamic parameters. IOP and BP were also measured four times, every 3 h from 9:00 a.m. to 6:00 p.m. over a single day or on two consecutive days ( Supplementary Fig. 2). That is, there were four measurement sessions (as specified above), and SS-OCTA, IOP, and BP were all taken for each session. To exclude any inter-measurement effects, all patients' examinations were taken in the order of systemic BP, IOP, and SS-OCTA. In addition, to minimize the residual effect of accommodation or pupillary dilation during the examination, each test interval was more than 5 min apart.
The IOPs of the recruited eyes were measured by Goldmann applanation tonometry in every session. All of the participants were in a sitting position and told to relax before measurement. IOP was measured three times for each session, and the median value was adopted. A single investigator (S.U.B.) performed all of the IOP measurements. The systemic BP including SBP and DBP was measured at the upper arm with an automated sphygmomanometer (OMRON HBP-1300, OMRON Healthcare, Co., Ltd China.) while the patient was asked to remain calm in a sitting position for at least 5 min.
According to the measured IOP and BP data for each session, the mean arterial pressure (MAP) and MOPP were analyzed. The MAP was calculated as follows: MAP = DBP + [1⁄3 × (SBP − DBP)] 35,36 . It was thus possible to calculate the MOPP at any specified time from the difference between the MAP and the IOP (substituted for venous pressure) as follows: MOPP = 2⁄3 × (MAP − IOP).

Intra-and
Inter-visit repeatability of RVD measurements. Repeatability of RVD was evaluated for the first 10 individuals of each group. The intra-visit repeatability of the peripapillary and macular RVDs was calculated from a subset of 10 healthy eyes and 10 low-teens NTG eyes, three sets of scans having been performed on a single visit. The ICCs, as obtained from three sets of scans performed on three separate visits, were used to evaluate test-retest variability 37 . The inter-visit reproducibility of the peripapillary and macular RVDs was analyzed from the same healthy subjects and low-teens NTG patients. Three OCT scans were performed within 30 min on 3 different days. The dataset measured for reliability was used independently from the dataset used for diurnal variation evaluation.
Statistical analysis. The sample size was determined by referring to existing studies 38,39 . The sample-size calculation was based on the assumption of a 10% difference in variation of RVD between the low-teens NTG and healthy groups with a 1-to-1 ratio of glaucomatous to healthy eyes. For a power of 90%, the required sample size was established to be 30 individuals per group.
Test-retest variability was evaluated by ICCs obtained with a two-way random model. The ICC values ranged from 0 to 1, a higher value indicating better reliability. ICCs < 0.40 were considered poor, those between 0.4 and 0.6 fair, those between 0.61 and 0.8 good, and those greater than 0.8, excellent 40 . Repeated-measures analysis of variance (ANOVA) was used to evaluate diurnal changes in RVD, IOP, BP, and MOPP in the two groups. Additionally, to investigate inter-group differences in diurnal changes of RVD, we ran a linear mixed model with adjusted ophthalmological and systemic factors followed by post hoc testing using the independent t-test. All other statistical analyses were performed using the Statistical Package for the Social Sciences version 22.0 for Windows (IBM Corp., Armonk, NY, USA). Statistical significance was defined as P < 0.05.

Data availability
Data supporting the findings of the current study are available from the corresponding author on reasonable request.